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Many extensions of the Standard Model include the possibility of light new particles, such as 
light Higgs bosons or dark matter candidates. These scenarios can be probed using the large 
datasets collected by B factories, complementing measurements performed at the LHC. This review 
summarizes recent searches for light New Physics conducted by the BABAR and Belle experiments. 
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I. INTRODUCTION 
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From supersymmetry to dark matter, many extensions of the Standard Model (SM) include the possibility of light 
^0 " New Physics. Thanks to their large luminosities, B factories offer an ideal environment to explore these theories. 
During the last decade, the BABAR Collaboration at PEP-II [l| and the Belle Collaboration at KEKB [2J, y] have 
respectively collected about 550 fb _1 and more than 1 ab _1 of data at several Y resonances, mostly the Y(AS) 
resonance (see Table fl]). These datasets have been exploited to explore many aspects of precision physics, including 
searches for light new particles. In the following, we review searches for light Higgs bosons, dark matter candidates, 
hidden sectors, sgoldstinos and Majorana neutrinos. 





BABAR 


Belle 


Total 


Y(5S) 




121 


121 


Y(4S) 


433 


711 


1144 


T(3S) 


30 


3 


33 


Y{2S) 


15 


25 


40 


Y(1S) 




6 


6 


Off-resonance 


54 


94 


138 



ov 

TABLE I: Integrated luminosities (fb 1 ) collected by the B factories at different center-of-mass energies. The off-resonance data 
^— ^ ■ were collected about 40 MeV below the T(4S) resonance at BABAR and at a similar offset for the T(4S) and Y(5S) resonances 
' in the case of Belle. 



II. SEARCH FOR LIGHT CP-ODD HIGGS BOSON IN Y DECAYS 

A light Higgs boson is predicted by several extensions of the Standard Model, such as the Next-to-Minimal Super- 
symmetric Standard Model (NMSSM). The NMSSM Higgs sector contains a total of seven states, three CP-even, two 
CP-odd, and two charged Higgs bosons. A CP-odd Higgs boson (A ) lighter than 2m?, can evade present experimental 
constraints making it accessible through radiative Y(nS) — > jA° decays 3. The corresponding branching fraction 
could be as large as a fewxlO -4 , well above the sensitivity of B factories [J, Q. 

The Higgs boson decay pattern depends on its mass and couplings, as well as the NMSSM particle spectrum. In the 
absence of light neutralinos, the A decays predominantly into a pair of muons below 2m T , while r + r~ and hadronic 
final states become significant above this threshold. The branching fraction A a — > x°X° may be dominant if the 
neutralino (x°) is the lightest stable particle with m x o < [7|. In this case, the neutralino is a natural dark 

matter candidate. 

BABAR has performed searches for a light CP-odd Higgs boson in a variety of decay channels. These measurements 
are discussed in the next paragraphs, and the results are summarized in Table [TTJ They place stringent constraints on 
light CP-odd Higgs models. 



A. Search for T(2S, 35) -> jA°, A 



The Y(2S, 35) — > 7^4°, A — > /i + /J~ candidates are reconstructed by combining a photon with a pair of oppositely- 
charged tracks. The energy of the photon in the Y center-of-mass (CM) frame is required to be greater than 0.5 GcV 
and one or both tracks must be identified as muons by particle identification algorithms. Events containing additional 
tracks and photons are rejected. The T(2, 35) candidates are then fit, constraining their CM energies to the total 
beam energy, and imposing a common vertex for the tracks. A series of unbinned likelihood fits to the dimuon mass 
distribution is performed to extract the signal. No evidence of A is observed; 90% confidence level (CL) limits on the 
branching fractions are established at the level of (0.26 — 8.3) x 10~ 6 for 0.212 < m^o < 9.3 GcV [8]. The limits as a 
function of the A° mass are shown in Fig.[TJ together with limits on the product B(A° — > /i + /i~)/y, where fy- denotes 
the effective coupling of the b quark to the Y meson Q . Slightly more stringent constraints have been recently 
derived by the BES-III Collaboration for some A mass hypotheses below 3GeV using J/ip —> 7/x + /i~ decays (Toj |. 
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FIG. 1: Left: 90% CL upper limit on the T(35) -> 7A , A -> 11+ V>~ (top) and Y{2S) -> jA°, A -> n+fj, 

(middle) as a function of the A mass derived by BABAR. The limits on the product B(A° — > H + l-i~)fr, are also shown 
(bottom). The J/tjj and ip(2S) regions (solid bands) are excluded from the search. Right: Product of branching fraction 
T(35) -> ~/A , A -» r+r" (top) and the corresponding 90% CL upper limit (bottom) as a function of the ^4 mass set by 
BABAR. The region corresponding to Xb.j(2P) — > yY(lS) transitions (shaded band) is excluded. 



B. Search for r(3S) -jA° , A -> r+r^ 

The two taus of the T(35) — > 7A , A — > r + r _ decays are identified through their leptonic decays, t + —5- e + v e v T 
and t + — > /i + ^j// T . The signal signature consists of exactly two oppositely-charged tracks, identified as muons or 
electrons, and at least one photon with an energy greater than 100 MeV in the CM frame. A set of eight kinematic 
and angular variables are used to further suppress the background, which arises mainly from radiative r production 
and two-photon processes. The signal yield is extracted as a function of m^o by a simultaneous fit to the photon 
energy distribution of the eej, /U/Z7 and e[i"/ samples. No excess is seen; 90% CL limits on the branching fraction are 
set at the level of (1.5 - 16) x 10~ 5 in the interval 4.03 < m A a < 10.1 GeV [Tlj |. as shown in FigQ] 



C. Search for Y(2S, 35) -> 7A , A -> hadrons 

The hadronic final states are selected from fully reconstructed the A° — > hadrons decays. The highest energy 
photon of the event is identified as the radiative photon from the Y(nS) decay; the A candidate is then constructed 
by adding the four-momenta of the remaining particles, constraining the A decays products to originate from the 
interaction point to improve the resolution. The signal yield is obtained by fitting the candidate mass spectrum in the 
range 0.3 — 7.0 GeV in steps of 1 MeV. The results are compatible with the null hypothesis; limits on the branching 
fraction are therefore set in the range (0.1 — 8) x 10~ 5 with 90% confidence level |12j |. 
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Mode Mass range ( GeV) BF upper limit (90% CL) 



Y(2S, 35) -> 7A , A -> 


0.21 < m A < 9.3 


(0.3- 


8.3) x 10~ 6 


r(35) ^ 7 a°,a ^t+t- 


4.0 < m A < 10.1 


(1.5- 


- 16) x 10 -5 


T(2S, 35) -» 7A , A -> hadrons 


0.3 < m A < 7.0 


(0.1 


- 8) x 10~ 5 


T(15)->7^ ,^°->XX 


m x < 4.5 GeV 


(0.5 - 


- 24) x 10~ 5 


T(15) -> 7 A , A -> invisible 


m A < 9.2 GeV 


(1.9- 


- 37) x 10~ 6 


r(35) -> 7 A , A -> invisible 


m A < 9.2 GeV 


(0.7- 


-31) x 10 -6 



TABLE II: Results of light Higgs boson searches performed by the BABAR Collaboration. 

D. Search for T(25) -> ir + ir~T(lS),r(lS) -> -yA°,A -> invisible 

This final state is characterized by a pair of low-momentum tracks, a single energetic photon, and large missing 
energy and momentum. Additional criteria on the photon and the extra neutral energy in the event are applied 
to further suppress contributions from electron bremsstrahlung, radiative hadronic T(1S) decays and two-photon 
processes in which particles escape undetected. The signal is extracted by a series of bidimensional unbinned likelihood 
fits to the dipion recoil mass and the missing mass squared for both two-body decays, T(1S) — > 7 A°, A — > invisible, 
and non-resonant three-body processes, T(1S) — > 7^4°, A — > XX . Values of m^o and m Xo are probed over < niA < 
9.2 GeV and < m Xo < 4.5 GeV, respectively. No significant signal is found; 90% CL limits B(Y(1S) -> -fA°,A° -> 
invisible) < (1.9 - 37) x 10~ 6 and B(Y(1S) -> 7^°,^° -> X °X°) < (1-9 - 37) x 1(T 6 are set 13]. These limits arc 
displayed in Fig. [2] as a function of the A and X ° masses. 

E. Search for Y(3S) -yA°,A° -> invisible 

The Y(3S) — > ~fA°, A — > invisible decays are also selected from events containing a single energetic photon with no 
additional activity in the detector. The background arises mainly from e + e~ — > 77, radiative Bhabha, and two-photon 
fusion events. The A yield is extracted by a series of unbinned likelihood fits to the photon energy distribution for 
< itiao < 7.8 GeV. No excess is seen, and limits on the branching fraction at the level of (0.7 — 31) x 10~ 6 are 
derived with 90% confidence level flij . 
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FIG. 2: 90% CL upper limits on the Y(1S) -> ^A°,A° -> invisible (left) and Y(1S) -> yA°,A° -> XX (right) branching 
fractions as a function of the A and x masses set by BABAR. 



III. SEARCH FOR DARK MATTER IN INVISIBLE T(15) DECAYS 

In a minimal model, a single dark matter particle ( X ) is added to the SM content, together with a new boson 
mediating SM-dark matter interactions [L5l - fl7l ]. A light mediator could be produced in bb annihilation and decay 
into a XX pah, contributing to the invisible width of Y mesons. In the SM, invisible Y(1S) decays proceed via 
the production of a vv pair with a branching fraction B(Y(1S) — > vv) ~ (1 x 10 -5 ) [l8j], well below the current 
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experimental sensitivity. The rate T(1S) — > XX is predicted to be larger by one or two orders of magnitude than that 
of T(1S) —> vv, assuming no flavor changing currents [19|. 

A search for dark matter in invisible T(1S) decays has been performed by BABAR using a sample of 122 x 10 6 T(3S) 
mesons [2(|. The T(1S) mesons are selected by reconstructing the T(3S) — > it + tt~T(1S) transitions. The dipion 
recoil mass peaks at the T(1S) for signal events, while backgrounds are broadly distributed. This strategy is common 
to several analyses, and provides a very clean T(1S) sample. 

The event topology consists of exactly two oppositely-charged tracks without any additional activity. The selection 
is performed using a multivariate classifier based on variables describing the pions, the neutral energy deposited in 
the calorimeters and the multiplicity of K\ candidates. 

The distribution of the resulting dipion recoil mass (Fig. [3]) shows a clear peak corresponding to T(1S) mesons 
on top of a non-resonant component. In addition to signal events, a background from T(1S) decays in which the 
decay products escape undetected, is also present. This component, kinematically indistinguishable from the signal, 
is evaluated using Monte Carlo simulations. 

The sum of signal and peaking background yields is first extracted by an extended maximum likelihood fit to the 
dipion recoil mass. After subtracting the peaking background, a signal yield of — 118 ± 105 ± 124 is measured, where 
the first uncertainty is statistical and the second systematic. No evidence for Y(1S) — > invisible decays is observed 
and a 90% confidence level Bayesian upper limit on its branching fraction is set at 3.0 x 10~ 4 using a prior flat in 
branching fraction. This result improves the best previous measurement [2ll | by nearly an order of magnitude, and 
sets stringent constraints on minimal light dark matter models. 




FIG. 3: The distribution of the dipion recoil mass (M rec ) for BABAR data, together with the result of the maximum likelihood 
fit (full line). The non-resonant background (dashed line) and the sum of the signal and the peaking background (solid filled) 
are also shown. 



IV. SEARCH FOR DARK MATTER AND HIDDEN SECTORS 



A new class of dark matter model has recently been proposed, following observation from satellite and ground- 
based experiments. These models introduce a new hidden sector with WIMP-like fermionic dark matter particles 
charged under a new Abelian gauge group [22H24| . The corresponding gauge boson, dubbed a hidden photon (A'), is 
constrained to have a mass at the GeV scale to explain the electron/positron excess observed by PAMELA [25[ and 
FERMI [26| , without a comparable anti-proton signal. The hidden photon couples to the SM photon through kinetic 
mixing with a mixing strength e, connecting the hidden sector to SM particles (27j . 

The Higgs mechanism generates the hidden boson masses .adding hidden Higgs bosons (h!) to the theory. A minimal 
model includes a single hidden photon and a Higgs boson [28|]. Additional gauge and Higgs bosons are considered in 
more complex variations [2!| [3(| ■ 
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A. Search for a hidden photon 



Hidden photons can be readily formed in e + e~ — > jA' interactions, and be reconstructed via their leptonic decays 
as resonances in the e + e~ — > -fl + l~ (I = e, u] spectrum. This signature is similar to that of light CP-odd Higgs 
production in e + e~ — > T(2S, 35 1 ) — > 7/z + /it _ and searches for this channel have therefore been reinterpreted as 
constraints on hidden photon production [3l[. The limits are shown in Fig. |4j together with bounds derived from 
measurements of <j> — > r]A' : A' —> e + e~ decays at KLOE (32[, and searches for direct production in fixed-target 
experiments (3ll . [33L A hidden photon could also contribute to the muon anomalous magnetic moment, and 

constraints from this measurement are also shown. Values of the mixing strength down to 10~ 3 — 10~ 2 are probed 
for 0.212 < m A > < 9.3 GeV. 

Other measurements could be reinterpreted as bounds on hidden photon production, such as searches for peaks 
in e + e~ — > 7t + t - events [ill ] or inclusive e + e~ — > 7 hadrons production [l2j . Invisible decays, occurring if hidden 
bosons decay to long-lived states, could also be detected as a mono-energetic photon peaks in e + e~ —> 7 + invisible 
events [i~4| . 




FIG. 4: Constraints on the mixing strength, e, as a function of the hidden photon mass derived from searches in Y(2S, 3S) 
decays at BABAR, (orange shading) and from other experiments (gray shading). The red line shows the value of the coupling 
required to explain the discrepancy between the calculated and measured anomalous magnetic moment of the muon [35[ . 



B. Search for hidden sector bosons 



Non-Abelian extensions of hidden sectors introduce additional hidden gauge bosons, generically denoted W , W", .... 
The phenomenology depends on the precise structure of the model, but heavy hidden bosons decay to lighter states if 
kinematically accessible, while the lightest bosons are metastable and decay to SM fermions via their mixing with the 
hidden photon (2t| [3(| ■ Non-Abelian hidden sectors could also accommodate inelastic dark matter (3(| if the mass 
spectrum contains nearly degenerated states. 

BABAR has performed a search for di-boson production in e + e~ — > A 1 * — > WW" , W —> l + l~ , W" — > l + l~ (I — e, u) 
events, where the bosons are reconstructed via their decays into lepton pairs (37j . The study has been performed in 
the context of inelastic dark matter models, searching for two bosons with similar masses. 

The signal signature consists of two narrow dileptonic resonances with similar masses carrying the full beam 
energy. This topology is quite unique; the only backgrounds arise from QED processes. The signal is extracted 
as a function of the average dileptonic mass in the range 0.24 — 5.3 GeV. No significant signal is found; limits on the 
e + e~ — > A'* — > WW" cross-section are derived. The results are translated into limits on the product a^e 2 , where 
a D = <?i>/47r and gu is the hidden sector gauge coupling constant. Values down to 10 -10 are probed, assuming nearly 
degenerate bosons. 
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C. Search for a hidden Higgs boson 



The Higgsstrahlung process, e + e~ — > A'h',h' — > A' A', offers another gateway to hidden sectors. This process 
is one of the few suppressed by only a single power of the mixing strength, and the background is expected to be 
almost negligible. The event topology is driven by the boson masses. While Higgs bosons heavier than two hidden 
photons decay promptly, they become metastable below this threshold and either produce displaced decays or escape 
undetected. 

A search for hidden Higgs boson in Higgsstrahlung production in the prompt decay regime has been conducted at 
BABAR, based on a data sample of 521 fb _1 [3a |. The measurement is performed in the range 0.8 < mw < 10.0 GeV 
and 0.25 < rriA' < 3.0 GeV, with the constraint mh> > 2mA'- The signal is either fully reconstructed into lepton or 
pion pairs (exclusive mode), or partially reconstructed (inclusive mode). The exclusive modes contain of six tracks, 
forming three hidden photon candidates with equal masses and a total invariant mass close to the e + e~ CM energy. 
The six pion final state has a significantly larger background than the other exclusive modes and is excluded from the 
search. Only two of the three hidden photons are reconstructed as dileptonic resonances for the inclusive modes. The 
remaining hidden photon, assigned to the recoiling system, must have a mass compatible with the Higgsstrahlung 
hypothesis. 

No significant signal is observed, and upper limits on the e + e~ — > A'h', b! — > A' A' cross section are set as a function 
of the hidden Higgs and hidden photon masses. These bounds are finally converted into limits on the product a^e 2 . 
The results are displayed in Fig. [SJ Values down to 10~ 10 — 10~ 8 are excluded for a substantial fraction of the 
parameter space probed. These limits are translated into constraints on the mixing strength in the range 10~ 4 — 10~ 3 , 
assuming ao = a — 1/137. 
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FIG. 5: The 90% confidence level upper limit derived by BABAR on the product a^f 2 as a function of the hidden photon mass 
(mA') for selected values of hidden Higgs boson masses (m h i). 



V. SEARCH FOR DIMUON DECAYS OF PSEUDOSCALAR SGOLDSTINOS 

The observation of three E + — > fi~ events with a dimuon invariant mass clustered around 214 MeV by the 
HyperCP Collaboration [39| triggered much discussion about the possibility of a new light state X produced in 
S + — > X, X — > p + pi~ decavs. Speculations about the nature of this state included a pseudoscalar sgoldstino [4(|, a 
hidden sector photon [H, Hl[ or a light Higgs boson [42| . Subsequent measurements in e + e _ collisions [1, |33| and pp 
interactions [43| excluded the light Higgs boson and hidden photon hypotheses. 

The Belle Collaboration performed a search for this state in B — s> K*°X,K*° — > K + it~ ,X — > \i + [i~ (B^.oj) and 
B — > p°X,p —> tt + tt~,X — > (B p ox) decays using a sample of 675 millions B°B° pairs [Hj]. The B mesons are 

reconstructed by combining two well-identified muons with a K + n~ or 7r + 7r~ pair. The signal candidates are identified 
using the beam-energy constrained mass M\, c = \jE'£ e — Pb an d the energy difference AE = Eb eam — Eb, where 
Ebeam denotes the beam energy and pb (Eb) the momentum (energy) of the B candidate in the e + e~ center-of-mass 
frame. Signal B K ,o x {B p o x ) events are selected in the region 5.27 < M bc < 5.29 GeV and -0.03 < AE < 0.04 GeV 
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(—0.04 < AE < 0.04 GeV). The corresponding dimuon mass distributions are shown in Fig. [6] No events are observed 
in the signal region. 

Lacking evidence for a pseudoscalar dimuon resonance at 214.3 MeV, 90% C.L. upper limits on the branching 
fraction B -> K*°X and B -> p°X are set at the level of 2.26 x 10~ 8 and 1.73 x 10~ 8 , respectively. They rule out 
models II and III of the sgoldstino interpretation of the HyperCP observation [45| . 
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FIG. 6: Distribution of the dimuon invariant mass for B — > (left) and B — > p°fi + fi (right) signal events for Belle 

data. The signal regions are shown as dashed bands in the inserts. 



VI. SEARCH FOR LEPTON NUMBER VIOLATION AND A MAJORANA NEUTRINO 

Lepton number is conserved in low-energy processes in the SM model, but can be violated in a number of New 
Physics scenarios, such as models containing Majorana neutrinos. In this case, the neutrino is its own antiparticle, 
and reactions changing the lepton number by two units become possible. The most sensitive searches have so far been 
based on neutrinoless nuclear double beta decays 0vf3(3 (46|, but the nuclear environment complicates the extraction 
of the neutrino mass scale. Processes involving meson decays, such as B — > h~l + l + (h — ir,K,D and I = e, fi), have 
been proposed as a possible alternative. The presence of a Majorana neutrino could mediate such a reaction, and 
would appear as an enhanced peak in the mass spectrum of the hadron and one of the leptons [U EH . 

BABAR has performed a search for the lepton number violating B + — > h~l + l + decays with h — it. K and I = e,/x, 
based on a sample of 471 ± 3 millions BB pairs 49]. The B meson candidates are reconstructed by combining a 
hadron with a pair of tracks identified as leptons from particle identification algorithms. The background is suppressed 
through boosted decision trees (BDTs) using variables describing the event shape and the B meson candidate. The 
number of B — > h~l + l~ events is ex tracted by a multidimensional likelihood fit of the BDTs response and the beam- 
energy constrained mass tties — \/s /4 — p* B 2 , where p* B is the momentum of the e + e~ CM frame. No evidence for 
such decays is observed, and limits on the corresponding branching fractions are set. 

A similar analysis has been conducted by Belle in B + — > D~e + e + , B + — > D~e + fi + and B + — > D~ /i + /i + decays, 
followed by a subsequent D~ — > K + tt~tt~ decay [50[. The analysis is based on a data sample of 772 million BB pairs 
collected at the T(4S) resonance. The B candidates arc identified using the energy difference, AE = Eb — Et, eam and 
m ES . The signal region is defined as 5.27 < m ES < 5.29 GeV and -0.055(-0.035) < AE < 0.035 GeV for the e+e+ 
and e + p + (p + p + ) final states. No signal events are observed and 90% CL limits on the branching fractions are set, 
assuming uniform three-body phase space distributions for B + D~l + l' + decays. 

The results are reported in Table UTTT and the limits on B + — ) h~l + l + decays as a function of the Majorana neutrino 
mass m v — mi+h- are also displayed in Fig. [7] The sensitivities of the B h~fi + fi + channels is similar to that of 
recent measurements from LHCb in the same channel (Fll . [52[ , and are an order of magnitude more stringent than 
the previous results for B + — > h~e + e + decays 



VII. CONCLUSION 



B factories have proved to be versatile machines, ideally suited to search for light New Physics over a wide range of 
processes. The next generation of flavor factories are expected to improve the sensitivity of these searches by one to 
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Mode BF UL (1(T 8 ) 


mode BF UL (10~ 6 ) 


B + -> 7r~e+e + 3.0 
B + ^ K~e + e + 2.3 
B + ^ir-n + n + 10.7 
B+ -> K-H+H+ 6.7 


B+ -> D"e + e + 2.6 
B+ -> D~e + fj, + 1.8 
B+ -> D-n + fi + 1.0 



TABLE III: Results of lepton number violation searches performed by the BABAR, (left column) and Belle Collaborations (right 
column) . The limits are given at 90% confidence level. 
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FIG. 7: 90% CL upper limits (UL) on the branching fraction B(B + — > h l + l + ) as a function of the h l + mass (m h -i+) derived 
by BABAR. 



two orders of magnitude, further constraining the parameter space of these theories. Many more results are to come 
in the near future, and will hopefully contribute to elucidate the nature of physics beyond the Standard Model. 
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